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Abstract 

We propose a new pattern of the neutrino mixing matrix which can be parametrized as the prod- 
uct of an arbitrary Hermitian matrix and the well-known tri-bimaximal mixing matrix. In this 
scenario, nontrivial values of the smallest neutrino mixing angle 9 13 and the CP-violating phases 
entirely arise from the non- unitary corrections. We present a complete set of series expansion for- 
mulas for neutrino oscillation probabilities both in vacuum and in matter of constant density. We 
do a numerical analysis to show the non-unitary effects on neutrino oscillations. The possibility 
of determining small non-unitary perturbations and CP-violating phases is discussed by measur- 
ing neutrino oscillation probabilities and constructing "deformed unitarity triangles". Some brief 
comments on the non-unitary neutrino mixing matrix in the type-II seesaw models are also given. 
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I. MOTIVATION 
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Recent solar [1], atmospheric [2|, reactor |3| and accelerator |4j neutrino experiments have 
convincingly verified the hypothesis of neutrino oscillations, which can naturally happen if 
neutrinos are slightly massive and lepton flavors are mixed. This discovery indicates the 
existence of new physics beyond the Standard Model (SM). Although it is easy to add either 
Dirac or Majorana neutrino mass terms to the SM, it is highly non-trivial to reveal the 
essential meaning behind such terms and find a natural and qualitative explanation of the 
smallness of neutrino masses. A complete theory of neutrino mass generation has been 
lacking and is eagerly desirable. 

A low-energy effective theory responsible for the generation of neutrino masses might give 
rise to slight violation of the unitarity of the neutrino mixing matrix. If three light neutrinos 
are mixed with other degrees of freedom (e.g., with the light sterile neutrinos [5j, the heavy 
Majorana neutrinos 6|, or the whole tower of Kaluza-Klein states in some models with extra 
dimensions {?]]), their 3x3 flavor mixing matrix V appearing in the SM charged-current 
interactions will in general be non-unitary. Therefore, the deviation of V from unitarity can 
serve as an indicator of new physics beyond the SM. 

A generic non- unitary neutrino mixing matrix V can be parametrized as V = H ■ V 1 , 
where if is a Hermitian matrix which can be written as 



(1) 



with a, b, c being real and (ij = 12, 13, 23) being complex, and V is a unitary matrix 
which is usually parametrized in terms of three mixing angles and one Dirac CP-violating 
phase [8] as 





fa 


K 12 


K 13 


H = 




b 


^23 






^23 


C / 



/ 



Vn 



CnC 



12^13 



Si oC 



S 12 C 23 C 12 S 23 S 13 e 



12^13 



C 12 C 23 ■ S 12 S 23' S 13 e 



s 13 e 



A 



V 



S 12 S 23 



C 12 C 23 S 13 e 



" C 12 S 23 



S 12 C 23 S 13 e 



S 23 C 13 



C 23 C 13 J 



(2) 



1 There is no unique way to parametrize a generic 3x3 flavor mixing matrix. For example, the non-unitary 
matrix V can also be expressed as V — A ■ V , where A is a lower (or upper) triangle matrix and V is 
a unitary matrix [9(. Although different parametrizations are mathematically equivalent, they may have 
their own advantages in describing different phenomena of neutrino physics. 
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Here we have omitted two Majorana CP-violating phases of V since they are irrelevant to 
neutrino oscillations. 

Given the lepton mixing matrix V, the constraints on the moduli of the elements of 



VV' have been deduced in Ref. 
oscillations and precision electroweak tests: 

/ 



10] by combining the experimental data on both neutrino 



It is then easy to find that a, b, c 



0.994 ±0.005 < 7.0 x 10~ 5 < 1.6 x 10~ 2 
<7.0xl0~ 5 0.995 ±0.005 < 1.02 x 10~ 2 
< 1.6 x to" 2 < 1.02 x HT 2 0.995 ±0.005 



(3) 
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< 3.5 x 10- 5 , |« 13 | < 8.0 x 10~ 3 , and \k 23 \ < 



5.1 x 10 3 should hold. We can further denote a, b, c as a = 1 + e a , b = 1 + e b and c = 1 + e c ; 
namely, 
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where e a , e b , e c are all real. Bounds on e a , e b and e c can also be obtained from Eq. (3): 
|ej < 5.5 x 10~ 3 , |e 6 | < 5.0 x lO" 3 and |e c | < 5.0 x 10~ 3 . 

The effects of non-unitarity of V on neutrino oscillations have been discussed in some liter- 



ature 



10. 
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121 ] . In particular, the authors of Ref. [12| have used the same parametriza- 
tion of V as given above and explored the effects of CP violation induced by those non- 
unitary complex parameters in neutrino oscillations. 

In this paper, we start from an intriguing point of view that the realistic neutrino mixing 
matrix V might result from a non-unitary correction to the well-known tri-bimaximal mixing 
pattern |13j]. The latter is compatible with current experimental data very well and can be 
derived from a number of flavor symmetries and their spontaneous or explicit breaking mech- 
anisms [ljj]. Instead of building a specific neutrino model to realize such a phenomenological 
conjecture, here we shall concentrate on the consequences of V on neutrino oscillations. 
In our new neutrino mixing scenario, V = H ■ V Q with V being given by 
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(5) 
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and H being an arbitrary Hermitian matrix as shown in Eq. (1) or (4). To be specific, the 
non-unitary neutrino mixing matrix V reads 

/ JL,._ J_„-. _, , 1 ^ 1 ; , J_,,-. ^. , \ 
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(6) 



It is clear that the parameters of H lead simultaneously to the unitarity violation and 
the deviation from V . The resulting smallest mixing angle 6 13 to be measured in reactor 
i> e — > z/ e oscillation experiments is attributed to the small parameters k l2 ari d ^13 • So are 
the CP-violating phases of V . 

The remaining parts of this paper are organized as follows. In sections II and III, we 
develop a complete set of series expansion formulas for neutrino oscillation probabilities 
both in vacuum and in matter of constant density, respectively, by taking account of the 
non- unitary mixing matrix V given in Eq. (6). In section IV, we discuss the possibility 
of determining some parameters of H by constructing the "deformed unitarity triangles". 
Section V is devoted to a short discussion about incorporating our parametrization of the 
non-unitary neutrino mixing matrix into a generic type-II seesaw model. Finally some 
conclusions are drawn in section VI. 



II. NEUTRINO OSCILLATIONS IN VACUUM 

Suppose that the non-unitary V in Eq. (6) describes the mixing between the neutrino 
fields in the mass basis and those in the flavor basis, 

V « = V ai v i > (7) 

where a = e,//,r and i = 1,2,3. The probability of neutrino oscillation u a — > v» (P a/3 ) 
can be derived in a similar way to that in the unitary case. The procedures of deriving 
the formulas for neutrino oscillation probabilities can be found in Ref. [loj]. Here we follow 
another way, in order to be concise. The derivation may be easily understood as follows. 

A typical neutrino oscillation process v a — > v„ can be divided into three parts [isj]: 1) u a 
being produced at the source through the charged- current interaction which can be denoted 
as W — > l a v a - Here, v a is a superposition of mass eigenstates za; 2) u i propagates from the 
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source to the detector; 3) Up (a superposition of ia) being catched by the detector through the 
charged- current interaction UpW — > lp. Therefore, the amplitude of the neutrino oscillation 
u a — > Up can be correspondingly divided into three parts: 

My a -> u p ) = A (W - U)Prop(^) A(u,W -> Z„) . (8) 



In the case of non-unitary neutrino mixing, it follows that A(W — > Z^za) = V*J a/ 

The factor l/ v /(V\/t) aa ensures that the total rate P(W -> Z" Q i/ Q ) = £\ -> = L 



Similarly, we have A{u$V — > Z^) = Vp i / y (W^)pp- The expression of Prop(z/j) is the same 
as that in the unitary case: Prop(zA) = exp(— imfL/2E l/ ). Finally, the amplitude of the 
neutrino oscillation u a — > Up is given by 

1 



(9) 



Then P a/3 , the probability of neutrino oscillation u a — > Up, is given by 
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where Ay = Am?.L/(4£7 v ) with Am?. = m\ - m 2 , and A", = BeWJV^Vft, = 
^[V^VpjV^jVp*^. We can further absorb the renormalization factor 1/ (VVt) aa (W^) pp 
into the redefinitions of A 1 ^ and J l Jp and rewrite the above equation as 



4i 21 ^ sin 2 A 21 - AA^p sin 2 A 31 - AAfp sin 2 A 32 



-2 J 2 J sin 2 A 21 - 2 sin 2 A 31 - 2 sin 2 A 32 , 



(11) 



where A% = A%/(VV% a (VV^)pp, J% = J%, <(VV% a (VV%p. 

The first term in Eq. (11) is the so-called "zero-distance" term. It means that at L = 
we have 



PUL = 0) 
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(12) 



Note that if a = (3, P a p{L = 0) = 1; namely, there are no "zero-distance" effects in 
the disappearance experiments. If a ^ /3, the oscillation probability P Q/3 (L = 0) is in 
general nonzero, that is the "zero-distance" effect. One can find that although it is nonzero, 
this term add only a tiny constant to the oscillation probability, and does not change the 
oscillatory behavior. In our scenario, we have P e „(L = 0) ~ 4|£ 12 | 2 , P eT {L = 0) ~ 4|£ 13 | 2 
and P^ T (L = 0) ~ 4|£ 23 | 2 . Another significant difference between the non- unitary and 
unitary cases is: if the mixing matrix V is non-unitary, there may exist 9 different Jarlskog 
invariants J % K corresponding to 3 different oscillation channels instead of a unique Jarlskog 
in the unitary case. 

If A 21 <C 1 is satisfied, we can expand Eq. (11) as 



P 
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-4 (i 2 i + iS) A 2 , - 4 [A% + Ag) sin 2 A 31 
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-AJ% (A 21 sin 2A 31 - 2A 21 sin 2 A 31 ) , (13) 

In our calculations, we find that although all the nine J l h are of 0{k i j), only J^ 2 + J^. is of 
0(k i j) while Jg 2 + and Jg 2 + are both of («?•). This observation means that the 
most sensitive way at short baseline neutrino oscillation experiments to detect CP violation 



is to measure the v — > z/ T channel. Such a point was also pointed out in Refs. [9_|, [12]. 

Here we present a complete set of formulas for neutrino oscillation probabilities P af3 to 
the second order in powers of k l2 , k 13 , k 23 , e a , e b , e c and A 21 . These formulas are good 
approximations for the A 31 -dominated oscillations, i.e., for neutrino oscillation experiments 
with relative short baselines and relatively high energies. 

Pee « l-2|/c 12 + « 13 | 2 sin 2 A3 1 -^Ai lJ (14) 
P w « 1- [l-4(Re[/i 23 ]) 2 ]sin 2 A 31 + ^(l + 2Re[«: 12 ])A 21 sin2A 31 

-^(2-3sin 2 A 31 )A 2 1 , (15) 
Prr « l-(l-4(Re[/i 23 ]) 2 )sin 2 A 31 + ^(l-2Re[«: 13 ])A 21 sin2A 31 

-^(2-3sin 2 A 31 )A 2 1 , (16) 



P e » » 4|/t 12 | 2 - (3|£ 12 | 2 - |K 13 | 2 + 2Re[£ 12 ^ 3 ]) sin^ + ^A^ 
2 

+-Re[£ 12 + £ 13 ]A 21 sin2A 31 + 2Im[K 12 K^ 3 ] sin2A 31 
4 

+ - (2Im[/c 12 ] - lra[k 12 + « 13 ] sin 2 A 31 ) A 21 , (17) 
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2 

— o R- e [^i2 + ^13)^21 sm 2A 31 — 2Im[/c 12 £* 3 ] sin 2 A 31 

~ (2Im[« 13 ] - Im[£ 12 + « 13 ] sin 2 A 31 ) A 21 , (18) 



Put » 4|/l 23 | 2 + [l-|/l 12 + /l 13 | 2 -4(|«: 23 | 2 + (Im[/l 23 ]) 2 )]sin 2 A 



31 

- (1 + Re[£ 12 - k 13 ]) A 21 sin 2A 31 + ^ (l - 3 sin 2 A 31 ) A 2 X 

+ [( 2 -e b - e c ) Im[£ 23 ] - lm[k 12 k* l3 \] sin2A 31 
4 

-3 (2Im[/t 23 ] - Im[£ 12 + « 13 + £ 23 ] sin 2 A 31 ) A 21 . (19) 

The first term in each of the above six equations is the "zero-distance" term. The last two 
terms in Eq. (17), (18) or (19) are the "CP-violating" terms. 

Suppose that the absolute values of those non-unitary parameters are around their upper 
bounds, i.e., \k 12 \ ~ 3.5 x 1(T 5 , \k 13 \ ~ 8.0 x 1(T 3 , \k 23 \ ~ 5.0 x 10~ 3 , |ej ~ 5.5 x 10~ 3 , 
\e b \ ~ 5.0 x 10~ 3 and \e c \ ~ 5.0 x 10~ 3 . We notice that P ee is only sensitive to \V e3 \ = 
\k 12 + k 13 \/y/2, therefore we are able to determine \k 13 \ through the detection of v e — > v e 
oscillation. By measuring the probability of v — > v oscillation it is possible to determine 
or constrain the value of Re[/c 23 ]. If the small difference between P TT and P can be well 
measured, then we are able to obtain the information on Re[£ 13 ]. Combined with the value 
of I &i 3 \, arg(/=c 13 ) can be determined. As for Im[«; 23 ], the most effective way is to probe the 
CP-violating terms in the v — > v T channel. 



III. NEUTRINO OSCILLATIONS IN MATTER 



When a neutrino beam passes through matter, only v e can interact with electrons in 
the medium via the charged-current interactions, while v e , and v T can all interact with 
electrons, protons and neutrons in the medium via the neutral-current interactions. The 
coherent forward scattering from the constituents of matter modifies the evolution behaviors 
of the neutrino beam. In the vacuum mass eigenbasis, the evolution equation of neutrinos 
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can be written as 



(20) 



We use tildes to denote the quantities in matter. For the propagation of neutrinos in matter 
of constant density, the Hamiltonian ii is given by 



U = E + V 1 AV* 



(21) 



where E = diag(_E 1; E 2 , E 3 ) is the energy matrix in the mass eigenbasis in vacuum, A 

I 
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G F n n (n e and n n are 



diag(V ( 

the electron and neutron densities, respectively). Here V is just the non-unitarity mixing 
matrix in Eq. (6). 

The Hermitian matrix 7i can be diagonalized by a unitary transformation 7i = UEW, 
where E = diag(i? 1 , E 2 , E 3 ) is the effective energy matrix in matter. The solution to Eq. 
(20) can be expressed as 

W m (L)) = Ue-^\u m (Q)), (22) 

where we have inserted L = t. From Eq. (22) we can work out the neutrino oscillation 
probabilities in matter: 
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(23) 



where X = VU*, A% = Re[X ai X pj X* aj X^ J% = lm[X ai X 0j X* .Xft and A,. = 
Comparing Eq. (23) with Eq. (11), we find that the matrix X, which is also non-unitary, 
can be regarded as the effective neutrino mixing matrix in matter. 

In Appendix A, we present the details of the approximate diagonalization of the Hamilto- 
nian 7i by using the perturbation theory. In the results to be presented below, those terms 
in proportion to (V cc — 2V NC ) will be neglected. The reason is simple: for ordinary earth 
matter, which is electrically neutral, we have n e ~ n n to a good degree of accuracy, and 
thus we can safely set V cc — 2V NC = \/2G F (n e — n n ) ^0. It is necessary to mention that 
the subsequent analytical approximations are not very good for a relative large L/E u or for 



S 



the Am 21 -dominated oscillation. In addition, we cannot directly obtain Eqs. (14) ~ (19) 

from Eqs. (27) ~ (32) by setting V cc , V NC — > 0. This is because the expansion of H in Eqs. 

(A4) ~ (A7) is improper if V cc = V NC = 0. 

The eigenvalues of 7i, E i (i = 1,2,3), are related to the effective neutrino masses in 

A 2 

matter by the relations E i E v + where E v is the energy of neutrinos 2 . Then the 
effective mass squared differences in matter are given by AA?- = 2E v (E i — Ej) which are 
shown in Eqs. (24), (25) and (26) to the second order in k 12 , k 13 , k 23 , e a , e b , e c and A 21 . 
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In Eqs. (27) to (32), we present the expansion forms of Eq. (22) for all six neutrino 
oscillation probabilities to the second order in k 12 , k 13 , k 23 , e a , e b , e c and A 21 in terms of 



Here we use the notation X i and AA?- to denote the effective neutrino masses and the mass squared 
differences in matter instead of m i and Am| , because we did not ordering E i according to their magnitude 
and the mass spectrum. After ordering E i , we have \ = rh i and AA? = Am?- . 
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In order to obtain the probabilities of anti-neutrino oscillations u a — > i/g, one needs to 
simultaneously change the signs of V^. c , and the terms underlined in the expressions of 

Pa/3- 

Different from the case in vacuum, the terms of 0(k i j) which can be strongly enhanced 
by large L appear not only in the expression of P^ T but also in those of P and P TT . If we 
omit the terms of O(k^), then we get very concise formulas for P ' P TT and P^ T : 

P m = l-sm 2 A 31 + 2^A 21 + 2V NC LRe[k 23 ]ym2A 31 , (33) 

P TT = l-sm 2 A 31 + 2t^A 21 + 2V NC LRe[k 23 }^sm2A 31 , (34) 

P MT = sin 2 A3 1 + 2QA 21 + 21/ 7VC LRe[K 2 3]^sin2A3 1 + 2Im[/€ 2 3]sin2A 31 . (35) 

We carry out a numerical analysis to show the difference between the corrections from 

the non-unitary parameter k 23 and the corrections from the nonzero 6 l3 to the neutrino 

oscillation probabilities. We compare between two special cases: Case I, we consider a 

unitary and nearly tri-bimaximal mixing matrix with nonzero 8 13 , as V given in Eq. (2) 

with 6 l2 = arcsin —= and 9 23 = 45°; Case II, we consider the non-unitary mixing matrix as 
V3 

shown in Eq. (6). The inputs of our numerical calculations are Am 21 = 8.0 x 10~ 5 eV 2 , 
Am 31 = ± 2.5 x 10~ 3 eV 2 (sign "+" for the normal hierarchy and "— " for the inverted 
hierarchy), the matter density p = 2.7 g/cm 3 . And in Case I we choose 6 l3 = 10°, in Case 
II we choose e a = -5.5 x 1(T 3 , e b = -5.0 x 10~ 3 , e c = -5.0 x 10~ 3 , k 12 = 3.5 x 10~ 5 • e'i, 
k 13 = 8.0 x 10" 3 • e*4 and \k 23 \ = 5.0 x 10~ 3 . Our numerical analysis is independent of our 



analytical results given above, but it confirms the results of our ana 



As for the analytical approximations for Case I, one may refer to 
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ytical approximations. 
17\ , in which there exists 



a complete set of series expansion for three-flavor neutrino oscillation probabilities in matter 
in terms of small 8 13 and a = Am^/Am^. 

Fig. 1 shows the effective mass differences in matter as functions of the neutrino beam 
energy E v in Case I and Case II for both the normal and the inverted hierarchies. We 
can clearly see that the mass squared differences are strongly magnified if E u is large (or 
equivalently if the matter density is large). An interesting point is that in Case II the effective 
mass difference Arh| 2 can reach zero at around E u — 12 GeV in the normal hierarchy case 
while in Case I the nonzero sin 8 l3 ensures the nonzero value of Am 32 . If we choose a nonzero 
8 l3 for Case II, we will get similar curves as those of Case I. We find that although the non- 
unitary parameter —= (k 12 + k 13 ) plays a very similar role as sin 8 l3 e~ lS in the expressions 
V2 

of neutrino oscillation probabilities in vacuum, it has very different effects from sin 8 l3 e~ lS 
in matter. This finding provides us with a possibility of distinguishing the nonzero 8 l3 from 
the non- unitary parameters in V. 

Taking account of Eqs. (27) to (29) in Ref. [l^], we can easily see that the Dirac phase 5 
dose not appear in the expressions of the eigenvalues of 7i. However, in the non-unitary case, 
all the effective mass squared differences contain the terms proportional to E v V NC Ke[k 23 ] 
which is relevant to arg(£ 23 ) for fixed |/c 23 |. Fig. 2 shows the effective mass differences in 
matter as functions of the Dirac phase 5 in Case I or the phase of k 23 in Case II for both 
mass hierarchies, where we have chosen E v = 50 GeV. We find that the correction from the 
term E u V NC ~Re[k 23 ] can be around 10~ 5 eV 2 in this situation. 

Fig. 3 tells us how the probabilities of neutrino oscillations v„ — * v and v„ — > v T in 

71 

matter are modified with the changing of the baseline L, where we have chosen 5 = — and 
arg(£ 23 ) = — . We can clearly see from the figure that the probability P^ (P M/ J can be 
largely enhanced (depressed) by a long baseline if the matter effect is taken into account. At 
the baseline L = 4000 km, P is about 10~ 3 . Fig. 4 shows the probability P and P as 
functions of the Dirac phase S in Case I or the phase of k 23 in Case II for both mass hierarchies 
with the baselines L = 1000 km and L = 4000 km. We find that if P,„ can be measured 

[IT 

to the level of 10~ 4 at 4000 km from the source and I be well measured, arg(£ 23 ) 

may convincingly be determined. Another point worthwhile to point out is that from Eqs. 
(33) to (35) we find that these three probabilities have approximate sign[ATO 31 ] - arg(£ 23 ) 
degeneracy, which means P m TT> ^(Am^ arg(re 23 )) w P m TT> ^(-Am^ -arg(£ 23 )). This 
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degeneracy is broken by the term -A 21 sin2A 31 , which increases with the increase of the 
baseline L. Fig. 4 shows that this degeneracy breaking can reach 10~ 3 at L — 4000 km if 
the energy of neutrinos is taken to be E v = 50 GeV. 

We admit that it is very difficult to measure the transition probabilities to the accuracy 
of 10~ 3 or even 10 -4 in the present or forthcoming neutrino oscillation experiments. Given 
the small effects of unitarity violation, however, our numerical results can at least serve to 
illustrate how sensitive an ambitious long-baseline neutrino experiment should be to this 
kind of new physics. It is worth remarking two positive aspects of searching for the non- 
unitarity of V in the v disappearance or v T appearance experiments. First, the signatures 
of the non-unitarity can be strongly enhanced by the matter effects, and thus a high energy 
and a very long baseline are essential to detect appreciable effects of the non-unitarity in the 
neutrino oscillation experiments. As for neutrinos of energy around 50 GeV, a baseline much 
longer than 4000 km may have much better sensitivity to the non-unitary parameters in the 
neutrino mixing matrix, in which case the varying of the terrestrial matter density need to 
be taken into account [18| . Second, the dependency of the non-unitary parameters on 
the energy of the neutrino beam is different from other neutrino mixing parameters, and 
thus measuring the energy dependency of the transition probabilities will help to identify 
the small effects of the non-unitarity. 



IV. CONSTRUCTING "DEFORMED UNITARITY TRIANGLES" 

In the unitary limit, the neutrino mixing matrix V which relates the neutrino mass 
eigenstates (u 1: u 2 , v 3 ) to the neutrino flavor eigenstates (u e , v , v T ) is unitary. The unitarity 
implies 





Ki^Ti + v e2 v; 2 + v e3 v; 3 = 


o. 


V 


v el v T \ + v e2 v; 2 + v e3 v; 3 = 


o, 




^iK* + v, 2 v; 2 + v M3 v; 3 = 


o, 


A 3 : 


v el v; 2 + v^v; 2 + v Tl v; 2 = 


o. 


A 2 : 


v el v; 3 + v Ml v; 3 + v Tl v; 3 = 


0. 


A x : 


v e2 v; 3 + v, 2 v; 3 + v r2 v; 3 = 


. 



(36) 
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In the complex plane, Eq. (36) corresponds to six unitarity triangles [19| denoted as A T , 
A , A e , A 3 , A 2 and A l respectively. The area of each unitarity triangle is \\J\, where J 
is the Jarlskog invariant measure of CP violation for the unitary MNS mixing matrix. If 
there is no CP violation (e.g., the tri-bimaximal mixing), the unitarity triangles shrink to 
segments. In other words, introducing the "unitarity triangles" provides a geometric way to 
describe CP violation (which can be determined by the area of each triangle) by measuring 
the CP-conserving quantities (the sides of the triangles). 

If the mixing matrix V is non-unitary, the orthogonal relations in Eq. (36) are in general 
not satisfied and the unitarity triangles in the complex plane are in general open. Note that 
every two vector sides can determine one triangle in the complex plane, of which twice the 
area corresponds to one of the nine Jarlskog invariants. To be explicit, the Jarlskog invariant 
JaB ec l ua l s to twice the area of the triangle determined by V ai V^ and V^Vg ■ (or V ai V*j and 
VpjVpj). Apparently, if all the six triangles are closed, these nine Jarlskog invariants are all 
equal (to J). 

As for our special scenario, the CP-violating effects are attributed to the phases of k 12 , k 13 
and k 23 . Thus we are able to determine those unitarity- violating parameters by constructing 
those "deformed unitarity triangles". In Appendix B, we give all the eighteen sides of the 
six deformed triangles. Table I shows the ratios of two sides of A r , A^ and A e to the first 
order of k 13 and k 23 . Here we omit the smallest parameter k 12 . 

One can find from Appendix B that for A T we have S 1 + S 2 + S 3 ~ 2k l2 « 0, therefore 
A T is almost closed. From Table I we find that in our special scenario, if the ratio IS^I/jS^I 
of A^ can be well measured, I 

^13 1 Ccl-D- be determined. Combined with the ratio |5 , 2 |/|5' 1 |, 

the phase of 

^13 CcLIl also be obtained. The results can be and should be checked by the 
measurements of l^l/l^l and IS^I/ISJ of A , which is a validation of our scenario. In 
addition, constructing the triangle A e can give us the correlation between |k 23 | and arg(£ 23 ) 
if k l3 is fixed. 



V. FURTHER DISCUSSIONS 



In this section we would like to present a short discussion about the well-accepted see- 
saw mechanism [6J|, of which the unitary violation of the mixing matrix V is a general 
consequence. We show here how to incorporate our parametrization into a generic type-II 
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TABLE I: Ratios of the absolute values of the sides of A T , A^ and A e , to the first order in k 13 and 
k 23 . The smallest parameter k l2 1S omitted. 



A, 



A, 



S2I 
s i\ 



v e2 v; 2 



VelVZl 



7,7*, 

el fil 



\So 



13 



V e zV* 3 



' el nl 



"13 I 



v e2 v* 2 



7,7*1 

el rl 

7,7% 

v e3 y t3 



|1- 2«13l 



7,7*1 

el rl 



^3 

15, 



v e3 v; 3 



7,7*1 

el rl 



"13 1 



Si] 
5x| 



^3^*3 



\1 4k 23 + k 13 \ 



-|1 — 4k 23 — k 13 | 



211 - 3k 



13 1 



seesaw model. In the type-II seesaw model 21[ which contains n right-handed neutrinos, 
the neutrino mass terms can be written as 



1 

2 



M L M D 
ME M R 




+ h.c. , 



(37) 



where z/£ = Cv^ with C being the charge conjugation matrix, likewise for iV R . Here, 
M L is a 3 x 3 matrix, M D is a 3 x n matrix and M R is a n x n matrix. The overall 
(n + 3) x (n + 3) neutrino mass matrix in £ mass , denoted as A4, can be diagonalized by the 
unitary transformation WM.U* = Ai; or explicitly, 




M L M D 




M v 

E 



(38) 



where M v = diag (m 1; m 2 , rn 3 ) and M N = diag (M l5 ■ • • , M n ) with and Mj being the light 
and heavy Majorana neutrino masses, respectively. Note that the submatrices 7, U, R and 
S are all non-unitary. Suppose the mass eigenstates and the flavor eigenstates of the charged 
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leptons are identical. In the basis of mass states, the standard charged- current interactions 
of neutrinos can be written as 



C = J- 



w: 



\N nJ 



w: 



+ h.c. . (39) 



One can find that the matrix R describes the strength of the charged-current interaction 
between the charged leptons and the heavy neutrinos. In addition, the deviation of V from 



20] 



unitary is also characterized by R, since the unitarity of U requires VV' + RR' = 1 

In fact the diagonalization of M. can be divided into two steps by decomposing U into 
the product of two unitary matrices W and V: 



M L M D 
ME M Ki 



M v 
Mv 



V* 



M„ 

o E 



(40) 



where W and V take the general forms 



W 



ft 




!• -1 























(41) 



The matrices ?7 L , B, C and C/^ are in general non-unitary, but V L and V R are unitary 
matrices. We can easily find that: V = U L V L , R = BV R . Let's count the degrees of freedom 
of these matrices. The (n + 3) x (n + 3) unitary matrix contains (n + 3) 2 degrees of freedom. 



Suppose freedom left for W is (n + 3) 2 — 3 2 



n 



2x3 the parameters in V L and are 



22J is to suppose that B is an 



all free, the degrees of xn. An ansatz made for >V in Ref. 
arbitrary 3 x n matrix which contains just 2 x 3n degrees of freedom and then parametrize 
W as 



W 



-fit Vl - £tfi 



(42) 



Comparing Eqs. (38), (40), (41) and (42) with the parametrization V = H ■ V shown in 
Section I, we can simply choose: V = V L , then we can find that H = U L = Vl - BB^ = 
Jl- BV R V R 'Bi = Vl - RR ] which contains the same degrees of freedom as the Hermitian 
matrix in Eq. (1). That is to say in the parametrization we have chosen, the unitarity- 
violating parameters in H can be expressed as the functions of RR\ where R is of the order 
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of M D /M R . In other words, the deviation from the unitarity which is described by e in Eq. 
(4) is of the order of (jf^j ■ 

In the canonical seesaw scenarios, the light neutrino masses are attributed to the leading- 
order contribution M v = -M D M^Ml (type-I seesaw) or M v = M L — M D M^M^ (type-II 
seesaw). In order to obtain the light neutrino mass scale m u ~ 0.1 eV, the mass scale of 
right-handed Majorana neutrinos is expected to be as high as m R ~ 10 14 GeV. In such a 
case, one finds M D /M R ~ 1CT 12 or equivalently e ~ 10~ 24 , which is too small to be detected. 
The possible ways out have recently been discussed (see, e.g., Refs. 23] and |24(). For 



instance, one may first impose certain flavor symmetries on the textures of M D , M R and M L 
to guarantee M D M R 1 M'£ = (type-I seesaw) or M L — M D M R l = (type-II seesaw) and 
then introduce slight perturbations to them so as to produce the tiny light neutrino masses. 
As a consequence, the mass scale of three light neutrinos in this approach is essentially 
decoupled from the magnitude of R. In the models proposed in Refs. [23j, |24|, for example, 
the right-handed Majorana neutrinos are assumed to be around the TeV scale such that 
M D /M R ~ O (10 _1 ) holds and the elements of e can be as large as ~ O (10 -2 ), just close to 
their upper bounds constrained by current experimental data. These kinds of seesaw models 
will be tested at the LHC or ILC [251 ]; and one of their low-energy consequences, which is 
just the unitarity violation of V under discussion, will also be tested in the long-baseline 
neutrino oscillation experiments. 



VI. SUMMARY 



As we are about to enter an era of high precision neutrino physics, a general and com- 
prehensive analysis of the non-unitary neutrino mixing matrix makes sense and will be 
useful for phenomenological explanations of future measurements and tests of type-I and 
type-II seesaw models. In this paper we have investigated a new pattern of the neutrino 
mixing matrix which is the product of an arbitrary Hermitian matrix and the well-known 
tri-bimaximal mixing matrix. Starting with this non-unitary mixing matrix, we have pre- 
sented a complete set of series expansion formulas for neutrino oscillation probabilities both 
in vacuum and in matter of constant density. We have carried out a numerical analysis to 
emphasize the importance of matter effects in the measurements of the non-unitary param- 
eters and in distinguishing their effects from the effects induced by small 6 l3 . We find that 



18 



measuring the probability of v — > v T or v — > v oscillation with large neutrino energy 
(e.g., ~ 50 GeV) and a relatively long baseline (e.g., several thousand km) is a viable way 
to detect those non-unitary parameters. We have also discussed the possibility of determin- 
ing the small non-unitary perturbations and the extra CP-violating phases by constructing 
the "deformed unitarity triangles" . Finally we have shown that our parametrization of the 
non-unitary mixing matrix can be naturally incorporated into the generic type-II seesaw 
model. 
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APPENDIX A: DIAGONALIZING THE HAMILTONIAN IN MATTER BY US- 
ING THE PERTURBATION THEORY 

In this appendix we use the perturbation theory to diagonalize the effect Hamiltonian 
7i = U EW in matter of constant density. In the series expansion of H, we regard k 12 , 
k 13 , k 23 , e a , e b , e c and a = Am^/Am^ as the small parameters of the same order, and 
perform the diagonalization to the second order of them. Now we are going to diagonalize 
the effective Hamiltonian H = E+ V T AV*. 

First we draw E l ■ 1 (which contributes only a pure phase e~ lElL to the oscillation am- 
plitude and nothing to the oscillation probabilities) out of 7i and rewrite 7i as 

fi = E' + V T AV* = -j- • diag (0, a, A<) + V T AV* , (Al) 

2E V 

where A = diag (V cc — V NC , —V NC) —V NC ). We find that it is easier to perform the diago- 
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nalization of H' = V HV<} = V E'Vj + H T AH*, with 



and 

H T AH* 



V EVj = 



2E„ 



l + e) T A(l + eY 





f° 





o\ 




/ 


1 


1 


-l\ 




1 








a 












2 





1 


1 


+ 3 




1 


1 


-1 






1° 


1 


l ) 




V 


-1 


-1 


V 





(A2) 



V, 



cc 



J 



( 



\ 



\ 



( 



+ 



e* + e* 









-V NC (26* + e* 2 ) -V NC -1, 



(A3) 



The term —V NC ■ 1 can be neglected for the same reason as omitting the term E 1 ■ 1. Then 
we will diagonalize H' with H' = WE'W^ using the perturbation theory. We can easily find 
that E = E' + E 1 — V cc and U = VqW. The matrix X, which describes the mixing in 
matter, can be expressed as X = VU* = HV Vq W = HW, which is also non-unitary. 
We write 



(A4) 



where 



n 



,(0) 



A ^31 

4£„ 



n 



l 

3 



2E„ 




1 
1 

/ 1 

1 





1 
1 

1 
1 



\ 



(l 

\ 

\ ( 

+ Vcc 



(A5) 



-1 -1 1 



ze„ ft 12 



ft 



12 



ft 



13 













-2V, 



JVC 





K 12 


K 13 


K 12 




^23 


\ K 13 


^23 





(A6) 



H' (2) = V, 



cc 



( e 



6„K 



€ a K l2 e a K l3 



\ 



(1^12 



ft 



121 



ft 12 ft 13 



e a K 13 K 12 K 1Z 



ft 



131 
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(el 



NC 



n2\ 



13 1 



(e n + e h ) k* + k*k 



a T H) "-12 T n '13 n '23 V c a ~~t~ e J K 13 + K 12 K 23 



\ 



e a e b) ^12 ^1.3^9.3 e f> + I ^19 I + |& 



"13^23 



121 



b 23l 



"-23 

12 



12^13 
2 



.(A7) 



^( e a + e c) ^13 + ^12^23 ( e 6 + e e) ^23 + ^12^13 e c + I ^13 I +1^231 y 



For the eigenvalues and the eigenvectors, we also write E' i = E'f^ + E'\ ^ + E'\ "* and 
v i = vf^ + vf^ + f (for i = 1, 2, 3). The unitary matrix W = (i^, v 2 , v 3 ). 

ii'^ can be easily diagonalized. The eigenvalues and the eigenvectors of H'^ are: 



M 0) 



V cc , E> 2 



~ (o) Am|i _ 



(A8) 



and 



-i 0) = (1.0.0)' r._r = (:).—.- — 



,(o) 



,(0) 



1 1 

1 '71' 71 



T 



Then, the first and the second order corrections to the eigenvalues are given by 



E> 



\h> (1) 



11 ^ e> {0) - E' m ' 



and the corrections to the eigenvectors are calculated by 



,(!) 



(0) 



5 - E 



= E 



H> {2) 



(A9) 



(A10) 
(All) 



(A12) 



E' 



E 



f) (E>r-Efy\ 



,(0) 



1 

2 



E 



{e>T-e> 



) J 



• v 



(o) 



(A13) 



where H>™ = vf H' {n) vf\ 



Inserting Eqs. (A8), (A9) into Eqs. (A10) and (A12), we obtain 



E>T 
E'T 



1 Am? 



21 



+ 2e a (V cc -V NC ) , 



3 2E V 
2 Am 2 

3"^ 1 -^c(^ + e c + 2Re[«: 2 3]) , 



(A14) 
(A15) 



21 



E'T = -V NC (e b + e c -2Re[k 23 }) ; 



(A16) 



and 



2V2A 



21 



2A 



3V CC L 







\ 



21 



1 V( (e b -e c + 2ImM (e 6 - c c - 2Im[£ 23 ]) 



2v/2A 



31 



3^ CC L 2v^A 31 

^ C (e b -e c + 2ImM) (e 6 - e c - 2Im[« 23 ]) 



31 



2V2A 



31 



(A17) 

where A^- = Amf-L/ (AE U ) with Am.. = m 2 — m| (for = 21,31,32), and the terms 
proportional to (V cc — 2V NC ) are omitted. 

We can further calculate X = HW to the first order in k 12 , k 13 , k 23 , e a , e b , e c and a: 

( n , _ 2v^A 21 



A 



l + e„ 



+ 



K 12 K 13 



K 12 + K 13 



\ 



2A 



V2 



V2 



3y cc L 

2j ^-^23 , ^vc( e 6- e c + 2Im ^2 3 ]) b + k 23 V NC (e b - e c - 2Im[£ 23 ]) 

~r K 12 7= h 



2A 



21 



+ k 



A 

c — k 



2V2A 



11 



V2 



2 V / 2A, 



31 



2.3 



13 



V2 



+ 



^vc( e 6- e c + 2Im [^23]) c + £ 23 V Arc (e b -e c -2Im[K 23 ]) 



2v/2A, 



31 



^2 



+ 



2V2A 31 
(A18) 



In this paper we do not order the eigenvalues of 7i according to their magnitude and 
the mass spectrum. This ordering does not change the oscillation probabilities. In this 
appendix we give the results of E'\ E'\ W<® and which are enough for calculating; 



~ (2) 

the probabilities P ee , P and P eT to the second order. E' \ and vf\ which are not shown 



,( 2 ) 



eel e/i 



here, only correct the terms of the second order in P , P^ T and P T 



APPENDIX B: SIDES OF SIX "DEFORMED UNITARITY TRIANGLES" 
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Sides of six "deformed unitarity triangles" to the first order in k 12 , k 13 , k 23 , e a , e b and e c . 
• A.: 

s, = Ki^Ti 



11 1 

-3-3 (e + e b) + g (5«i2 - «ia + 2 ^2 3 ) 



s 2 = v e2 v; 2 



v e3 v; 3 



11 1 

3 + 3 (e a + ej + 3 (2k 12 - k 13 - k* 23 ) , 

2 \ K 12 + K 13> ) 



(Bl) 
(B2) 
(B3) 



with S ± + S 2 + S 3 ~ 2k 12 . 



22 



5 1 = V el V T \ » I + I ( £a + e J + I (5k 13 - £ 12 - 2^ 3 ) , (B4) 

5 2 = V e2 V T * 2 » -1 - i (e o + e c ) + i (2« 13 - k 12 + ^ 3 ) , (B5) 

5 3 = V e3 V T * 3 » 1 (« u + £ 13 ) , (B6) 



with + + 5*3 ~ 2/c 13 . 



e" 



= V^V T \ « -^-^(e 6 + e c ) + ^(« 23 + ^2-«i3) , (B7) 

^2 = ^ 2 K*2 « - £ (C 6 + O + l (2^23 - «U + ««) » (B8) 

= V% « £ + £ (C 6 + O + >%3 , (B9) 



with S ± + S 2 + S 3 ~ 2k 23 . 



• A 3 : 



^ = V el V: 2 » ^(l + 2eJ-^=(Re[« 12 -« 13 ]-Mm[« 12 -« 13 ]) , (BIO) 



^2 = V^lft « -^(l + 2eJ + ^=(Re[« 12 + 2« 23 ]-3nm[« 12 ]) , (Bll) 



5 3 = K1K2 « -^(l + 2eJ-^=(Re[« 13 -2« 23 ]-Mm[«; 13 ]) , (B12) 
with S , 1 +5 2 +5 3 « — — (2e a - e fe - e c ) + — (Re[2/t 23 + « 12 - « 13 ] - 3ilm[« 12 - « 13 ]). 



A 2 : 



^ = « ^= («« + *Is) , (B13) 



^2 = ^1^3 « -^(l + 2e b ) + -^=(^2 + ^23]) , (B14) 



S 3 = V Tl V T % » -L(l + 2e c ) + i=(^ 3 + ilm[«23]) , (B15) 

1 2 
with 5 X + S 2 + S 3 « ^= (e 6 - e c ) - (£* 2 + «J 3 + iIm[K 23 ]). 



^ = K2K3 « ^= (^2 + «*i3) , (B16) 
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V V* 

V t2 V t3 



Vq 
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with S'j + $2 + S3 
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FIG. 1: The effective mass-squared differences in matter as functions of the neutrino beam energy 
E v in Case I (the unitary case, represented by dashed lines) and Case II (the non-unitary case, 
represented by solid lines) for both the normal (the first plot) and the inverted (the second plot) 
hierarchies. 
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FIG. 2: The effective mass-squared differences in matter as functions of the Dirac phase 5 in Case 
I (the unitary case, represented by dashed lines) or the phase of k 2 ^ in Case II (the non-unitary 
case, represented by solid lines) for both mass hierarchies, where we choose E v = 50 GeV. 
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FIG. 3: The probabilities of neutrino oscillation v — > v (the first plot) and f„ — > v T (the second 
plot) in matter as functions of the baseline L in Case I (the unitary case, represented by dashed 
lines) and Case II (the non-unitary case, represented by solid lines), in the normal hierarchy case. 
Here we choose E v = 50 GeV. The dotted lines in the figure show the corresponding probabilities 
if the neutrino mixing is the exact tri-bimaximal mixing. 
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FIG. 4: The probabilities of neutrino oscillation 1/ — > 1/ (a) and i/„ — ► i/ r (b) in matter as functions 
of the Dirac phase 6 in Case I (the unitary case, the dashed lines for the normal hierarchy, the 
dotted lines for the inverted hierarchy) or the phase of k 2 ^ in Case II (the non-unitary case, the 
solid lines for the normal hierarchy, the dot-and-dash line for the inverted hierarchy), where we 



choose E„ = 50 GeV. 
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